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Italian Text Pursuant to Art 14.2 



"ARCHITETTURA SCALABILE DI CALCOLO PARALLELO BELLA 
TRASFORMATA COSENO DISCRETA DI BLOCCHI DI PIXEL DI 
DIMENSIONI DIVERSE E COMPRESSIONS PER CODIFICA 
5 FRATTALE" 

Campo pi applicazione 

L'invenzione conceme in generale i sistemi di elaborazione di immagini, da 
registrare e/o trasmettere e piu in particolare a sistemi di compressione e codifica 
di immagine basati sul calcolo della trasformata coseno discreta (DCT) di blocchi 
10 di pixels di immagine. 

L'invenzione e particolarmente utile in codification video secondo lo standard 
MPEG2 pur essendo applicabile anche ad altri sistemi. 

Background deli/invenzione 

II calcolo della trasformata coseno discreta (DCT) di una matrice di pixels di 
15 un' immagine rappresenta un passo fondamentale di elaborazione dei dati di 
immagine sui risultati della trasformazione DCT e quindi operata una divisione 
per una matrice detta di quantizzazione alio scopo di ridurre piu o meno 
drasticamente Fampiezza dei coefficienti DCT, come presupposto alia 
compressione dei dati che awiene in fase di codifica, secondo un certo protocolio 
20 di trasferimento dei dati video da trasmettere o memorizzare. 

II calcolo di trasformazione coseno discreta e tipicamente eseguito su blocchi o 
matrici di pixels in cui e suddivisa un' immagine. 

Esigenze di velocita dei sistemi di elaborazione di immagini per la loro 
registrazione e/o trasmissione, impongono Fimpiego di architetture hardware per 
25 velocizzare i vari passi del processo di elaborazione tra i quali, eminentemente, ii 
passo di trasformazione coseno discreta per blocchi dei pixels di immagine. 
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L'implementazione hardware impone peraltro la predefinizione di alcuni 
parametri fondamentali, eminentemente della dimensione della matrice o blocco 
di pixels in cui suddividere un'immagine per le esigenze di elaborazione. 

Tale predefinizione pud rappresentare un pesante vincolo che limita le possibility 
5 di ottimizzazione di un sistema di elaborazione, ad esempio di un codificatore 
MPEG2, o la sua adattabilita a condizioni di impiego molto diverse tra loro in 
termini di requisiti di prestazione. 

E inoltre evidente Fenorme vantaggio sia in termini di riduzione dei costi che la 
disponibilita di un sistema integrato di elaborazione dei dati che possa essere 
10 programmato a calcolare in parallelo la DCT su piu blocchi di pixels di 
dimensione selezionabile tra un certo numero di dimensioni disponibili. 

SCOPI E SOMMARIO DELL'lNVENZIONE 

E evidente il bisogno e/o Putilita di disporre di un metodo ed un dispositivo 
hardware di calcolo della trasformata coseno discreta (DCT) in parallelo su piu 
15 blocchi di pixels, la cui architettura consenta una scalabilita della dimensione del 
blocco di pixels, ad esempio tale da consentire il calcolo della trasformata coseno 
discreta (DCT) per blocchi di 8x8 pixels, o di quattro blocchi di 4x4 pixels in 
parallelo, o di sedici blocchi di 2x2 pixels in parallelo, operando una selezione di 
dimensione del blocco (size). 

20 La scalabilita della dimensione del blocco di pixels e la possibility di eseguire il 
calcolo della trasformata coseno discreta (DCT) in parallelo su piu blocchi di 
dimensione frazionaria rispetto alia dimensione massima di blocco mediante una 
struttura hardware, e inoltre strumentale alFimplementazione di schemi ed 
algoritmi di compressione di immagini "ibridi", ed altamente efficienti. 

25 Ad esempio, in virtu della scalabilita della dimensione di blocco e della capacita 
di calcolo parallelo della DCT su piu blocchi secondo la presente invenzione e 
implementabile su una codifica frattale applicata nel dominio della trasformata 
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coseno discrete (DCT) anziche nel dominio dello spazio dei dati di immagini, 
come e di consuetudine. 



Rappresenta pertanto un aspetto importante deirinvenzione anche un nuovo 
metodo di compressione e codifica dei dati di un' immagine reso praticamente 
possibile dairimpiego di una struttura hardware di calcolo della trasformata 
coseno discreta (DCT) su blocchi di dimensione scalabile deirinvenzione e 
consistente essenzialmente nel 

• suddividere un'immagine definendo due distinti tipi di blocchi di suddivisione: un 
primo tipo, di dimensione N/i*N/i, denominati blocchi range essenzialmente non 
sovrapponibili gli uni agli altri ed un secondo tipo, di dimensione N*N 9 denominati 
blocchi domain traslabili per intervalli di N/i pixels e sovrapponibili gli uni agli altri 
(cioe traslandp sulla immagine originate la finestra che individua un blocco domain di 
un intervallo equivalente alia dimensione orizzontale e/o verticale di un blocco 
range); 

15 • calcolare parallelamente la trasformata coseno discreta (DCT) dei 2 1 blocchi 
range e di un relativo blocco domain; 

• classificare i blocchi range trasformati in funzione della loro relativa 
complessita rappresentata dalla somma dei valori della terna di coefficienti AC; 

applicare la trasformata frattale del dominio della DCT ai dati dei blocchi range con 
classificazione di complessita superiore ad una soglia prestabilita e memorizzare solo 
il coefficiente DC dei blocchi range con complessita inferiore a detta soglia, 
identificando un relativo blocco domain di appartenenza del blocco range in 
trasformazione tale da produrre la migliore approssimazione frattale detlo stesso 
blocco range e di cui e calcolata la trasformata coseno discreta; 

25 • calcolare un' immagine differenza tra ciascun blocco range e la sua 
approssimazione frattale; 

• quantizzare detta immagine differenza nel dominio della DCT impiegando una 
tabella di quantizzazione predisposta in funzione delle caratteristiche della vista 
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umana; 



• codificare Timmagine differenza quantizzata mediante un metodo di codifica 
basato sulla probability dei coefficienti di quantizzazione; 

• memorizzare o trasmettere il codice di codifica per ciascun blocco range 



DCT, di coefficienti AC indicano la complessita del blocco. Considerando il caso 
in cui N=8, i quattro coefficienti in alto a sinistra del blocco sono: 

10 • il coefficiente DC che occupa la posizione 00; 

• i tre coefficienti AC che occupano rispettivamente le posizioni 01, 10 e 11. 

I coefficienti AC sono usati per decidere a quale classe appartenga un certo blocco 
range: se la somma dei loro valori assoluti e inferiore a una determinata soglia T, 
il blocco range in questione viene definito "a bassa attivita"; se invece questa 
15 somma e uguale o maggiore di T, il blocco range viene definito "ad alta attivita". 

Per un blocco range a bassa attivita, i coefficienti AC sono piccoli e pertanto 
possono essere omessi senza influenzare molto la fedelta: in questo caso il blocco 
pud essere approssimato memorizzandone soltanto il coefficiente DC. 

Per un blocco range ad alta attivita, il proseguimento di codifica frattale 
20 dell'invenzione consiste nel ricercare due trasformate lineari idonee (per esempio 
rotazioni, ribaltamenti, ecc), se r, <f> , e un blocco domain, del quale la DCT e 
denotata con Fd(u,v), che soddisfino, almeno approssimativamente, la seguente 
equazione: 



5 



compresso nel dominio della DCT e il coefficiente DC di ciascun blocco range 



non compresso. 



Indicando con F R (u,v), la DCT di un generico blocco range, nel dominio della 




a) 
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Avendo identificato il blocco domain piu simile o omologo al blocco range in 
elaborazione, se ne memorizzano i parametri F D (u, v), r,<fr . 

Si calcola quindi l'immagine differenza tra il blocco range e la sua posizione 
frattale 

5 D(u,v) = F r (u,v)-t o <f> (F D (u,v)) (2) 

Quantizzando Fimmagine differenza D(u, v) si ottiene 




D Q (u, v) = INTEG[D(u, v) / Q(u, v)] (3) 

dove: 

• Dq(u,v) e, nel dominio della DCT, l'immagine differenza quantizzata; 

10 • Q(u, v) e una tabelle di quantizzazione progettata tenendo conto delle 
caratteristiche della vista umana; 

• INTEG e una funzione che approssima il suo argomento airintero piu vicino. 

Dopo la quantizzazione, la maggior parte dei coefficienti di Dq(u, v) sono nulli. E 
facile progettare una codifica, ad esempio la codifica di Huf&nan, secondo le 
probability dei coefficienti. Alia fine si memorizza il codice da memorizzare o 
trasmettere. Quando ogni blocco range e stato codificato, la procedura di 
compressione e terminata. 

L'oggetto deU'invenzione e piu chiaramente definito nelle rivendicazioni annesse. 

Breve descrizione dei disegni 

20 I diversi aspetti implementativi deirarchitettura scalabile del calcolo della 
trasformata coseno discreta delPinvenzione nonche di un efficiente metodo di 
compressione e codifica frattale da essa permesso, saranno piu facilmente 
compresi attraverso la seguente descrizione dettagliata di una forma di 
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realizzazione dell'architettura delPinvenzione e delle diverse modalita di 
funzionamento secondo una selezione di dimensione (size) del blocco di pixels di 
suddivisione deirimmagine e facendo riferimento ai disegni allegati, nei quali: 

la Figura 1 e uno schema a blocchi di un codificatore attuante una compressione 
5 "ibrida" basata su codifica frattale e DCT, secondo la presente invenzione; 

la Figura 2 e il grafo di flusso di calcolo parallelo della DCT di sedici blocchi di 
dimensione 2x2; 

la Figura 3 illustra Farchitettura per il calcolo di sedici DCT 2x2; 

la Figura 4 mostra I'ordinamento dei dati in ingresso per il calcolo di sedici DCT 
10 2x2; 

la Figura 5 mostra la fase PROCESS del procedimento di calcolo di sedici DCT 
2x2; 

la Figura 6 illustra Tarchitettura per il calcolo di quattro DCT 4x4; 

la Figura 7 mostra I'ordinamento dei dati in ingresso per il calcolo di quattro 
15 DCT 4x4; 

la Figura 8 mostra la fase PROCESS del procedimento di calcolo di quattro DCT 
4x4; 

la Figura 9 illustra Tarchitettura per il calcolo di una DCT 8x8; 

la Figura 10 mostra I'ordinamento dei dati in ingresso per il calcolo di una DCT 
20 8x8; 

la Figura 11 mostra la fase PROCESS del procedimento di calcolo di una DCT 
8x8; 

la Figura 12 mostra rarchitettura scalabile per il calcolo di una DCT 8x8 o di 
quattro DCT 4x4 in parallelo o di sedici DCT 2x2 in parallelo; 
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la Figura 13 mostra la fase INPUT per rarchitettura scalabile delPinvenzione; 



la Figura 14 mostra la fase PROCESS per rarchitettura scalabile dell'invenzione; 

la Figura 15 mostra i blocchi che implementano la fase PROCESS per 
rarchitettura scalabile delPinvenzione; 

5 la Figura 16 e uno schema di dettaglio del blocco QA; 

la Figura 17 e uno schema di dettaglio del blocco QB; 




la Figura 18 e uno schema di dettaglio del blocco QC; 
la Figura 19 e uno schema di dettaglio del blocco QD; 
la Figura 20 e uno schema di dettaglio del blocco QE; 
10 la Figura 21 e uno schema di dettaglio del blocco QF; 
la Figura 22 e uno schema di dettaglio del blocco QG; 

la Figura 23 illustra la fase ORDER per rarchitettura scalabile dell' invenzione; 

la Figura 24 illustra la fase OUTPUT per l'architettura scalabile dell' invenzione. 

Pur facendo riferimento in alcuni degli schemi illustrati nelle figure ad 
un'applicazione particolarmente significativa ed efficace deirarchitettura di 
calcolo parallelo su blocchi di dimensione scalabile della DCT che comprende una 
fase di compressione per codifica frattale dei dati di immagine, resta inteso che il 
metodo e rarchitettura di calcolo parallelo della trasformata coseno discreta 
(DCT) di una matrice bidimensionale di dati di ingresso per blocchi di dimensione 
scalabile, forniscono un'eccezionale liberta di implementazione di algoritmi di 
compressione particolarmente efficienti ed efficaci awalendosi della scalabilita e 
della capacita di calcolo parallelo della DCT. 
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Descrizione pi alcune forme pi realizzazione 



Le fasi di partizionamento e di calcolo della trasformata coseno discreta di una 
matrice bidimensionale di dati di ingresso verra ora descritta distintamente per 
ciascuna dimensione di blocco range, secondo un esempio di realizzazione, a 
5 partire dalla dimensione piu piccola di 2x2 dei blocchi di dimensione frazionaria 
per i quali il calcolo della DCT e eseguito in parallelo, fino alia dimensione 
massima di 8x8 del blocco piu grande. 

A questa descrizione seguira quindi la descrizione di un'architettura scalabile 
secondo necessita cambiando il valore della variabile globale size. 

10 II procedimento di calcolo parallelo di DCT delPinvenzione pud essere scomposto 
in diverse fasi: 

• Fase INPUT 

• Fase PROCESS 

• Fase ORDER 
15 • Fase OUTPUT 

che saranno descritte in successione per ciascun caso considerato. 

L'operazione di trasformazione DCT pud essere definita come segue. 

Per una matrice dei dati in ingresso Xn*n ~ i x i,j\® , la matrice di 



20 



uscita yw*N = [y m) „ \$ <m,n< N -I ,k definita da: 



2 . w xV V ( (2i + l)m \ f {2j + l)n \ 

.» = Jf 6 ( m M w ) Z E X U C ° V 2N ^ J C ° V 2N *) ' ( } 



dove 



1 

per n = 0 
\per\ <n <N-\ 
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Per convenienza si assume che N=2\ dove i e' un intero e / > 1. Si rimuovano 
e(m),e(n), e il valore di normalizzazione 2/N nell'Eq. (4), poiche e possibile 
reintrodurli in un passo successive Pertanto d'ora in poi si utilizzera la seguente 
versione semplificata delPEq. (4): 



V V f(2/ + l>7 1 \ {2j + \)m 1 



(5") 



Calcolo parallelo di sedici DCT 2x2 



Per N = 2 l'Eq.( 5) diventa: 



r(2z + l)n 1 r(2y + l>7 1 
/=oy=o L * J L * J 



(6) 



II grafo di flusso per una DCT 2x2 e mostrato in figura 2, nella quale A=B =C=J e 
10 i dati in ingresso e in uscita sono i pixel nelle posizioni (0,0), (0,1), (1,0), (1,1). 

Si consideri ora come calcolare in parallelo sedici DCT 2x2 in cui si suddivide un 
blocco 8x8. 



II procedimento e suddiviso in varie fasi, la cui veduta d'insieme e in figura 3. 

In questa figura vengono evidenziate le trasformazioni eseguite sul blocchetto 2x2 
15 costituito da (0,6), (0,7), (1,6), (1,7). 

I pixel che costituiscono il blocco in ingresso vengono ordinati nella fase INPUT 
ed elaborati nella fase PROCESS in modo da ottenere i coefficienti delle sedici 2- 
D DCT su 4 campioni. Per esempio, la 2-D DCT del blocchetto (0,1) costituito da 

{l[0],m[0],n[0],o[0]} e {a[0],b[0],c[0] t d[0]}. 

20 I coefficienti delle 2-D DCT vengono riordinati nella fase ORDER in otto vettori 
a otto componenti. Per esempio i coefficienti {a[0],b[0],c[0],d[0]} vanno a 
costituire il vettore F. 
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I vettori cosi ottenuti vengono inviati alia fase OUTPUT per ottenere i coefficienti 
della 

DCT 2x2, costituendo il blocco di uscita. 
Fase INPUT 

5 I pixel di ciascun blocchetto (i j), con 0 < / < 1 e 0 < j < 3 , vengono ordinati per 
costituire i vettori a 8 componenti l f m,n,o nel modo seguente: 

• i pixel che occupano nel blocchetto la posizione (0,0) costituiscono il vettore /; 

• i pixel che occupano nel blocchetto la posizione (0,1) costituiscono il vettore 
m; 

10 • i pixel che occupano nel blocchetto la posizione (1,0) costituiscono il vettore 
n\ 

• i pixel che occupano nel blocchetto la posizione (1,1) costituiscono il vettore 
o\ 

Analogamente i pixel di ciascun blocchetto (ij) 9 con 2 < / < 3 e 0 < j < 3 , vengono 
15 ordinati per costituire i vettori a 8 componenti p,q,r,s nel modo seguente: 

• i pixel che occupano nel blocchetto la posizione (0,0) costituiscono il vettore 

• i pixel che occupano nel blocchetto la posizione (0,1) costituiscono il vettore 

20 • i pixel che occupano nel blocchetto la posizione (1,0) costituiscono il vettore 
r; 

• i pixel che occupano nel blocchetto la posizione (1,1) costituiscono il vettore 

Questo ordinamento e mostrato in dettaglio in figura 4. 

25 Si osservi, per esempio, che i pixel del blocchetto (0,3) vanno a costituire la 
componente 3 dei vettori l,m,n,o. 
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Fase PROCESS 



La fase PROCESS consiste nel calcolare in parallelo le sedici 2-D DCT 
elaborando i vettori a otto componenti /,m,...,s come mostrato in figura 5. Si 
osservi, per esempio, che i coefficienti della 2-D DCT del blocchetto (0,3) vanno 
a costituire la componente 3 dei vettori a,b,c,ddi Figura 3. 

Fase ORDER 

La fase ORDER consiste nelPordinamento 

delle sequenze di uscita delle otto 2-D DCT in otto vettori l',m 9 , ...,5' cosi definiti: 



10 
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Si osservi, per esempio, che i coefficienti della 2-D DCT del blocchetto (0,3) 
vanno a costituire le componenti 4,5,6,7 del vettore m\ 
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Fase OUTPUT 



Questa fase consiste nel riordinamento dei dati in uscita: a partire dai vettori di 8 
componenti a,b,„,,h si costruisce il vettore di 64 componenti 

y cosi definito: 



y[5A\ y[55] 



y[63\ 



r/[o] 


m 


/[4] 


/[5] 


m[0] 


-[1] 


m[4) 




1[2] 
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/[6] 


/[7] 


m[2] 


43] 


m[6] 


«[7] 


n[0] 


n[l] 
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»[5] 
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o[l] 
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o[5] 


42] 


»[3] 
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HH 
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,[4] 
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K3] 


r[6] 
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,[2] 


s[3] 


,[6] 


,[7]_ 



(7) 



Calcolo parallelo di quattro DCT 4x4 

FerN = 4 V Eq. (5) diventa: 



3 3 

>v* = ZZ*/j cos | 

/=0 /=0 



(2i + 1>i 
S 



^■Jcos 



(8) 



Ponendo: 







>o~ 




"/o" 
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^16*1 = 


^1 
J>2 
_^3. 


• F y6*\ - 


A 

h 



dove: 



yi=[y.,o s yi,i 9 yu»yu]' 

Wr to = DCT U,i to\ to = ^fe,/ to 
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for t„ = scr(k<U. (/v to - ^({*uU 



{^3,/? =0 = Ko'*3,l>*0,2>*2,3}> 
fe,/?=0 = ^2,0» X 0,1»*3,2» X U ' 
{ B \,i }/=0 = ' ' ^L 2 ' X °. 3 



si puo dimostrare che: 



Y \6*\ = (^4)16^16*1. 



(9) 



dove 



(£4)16 = 



("0)4 
("1)4 
("2)4 
("3)4 



(#0)4 

(#3)4 
-(^2)4 

-(^1)4 



(H 0 ) 4 (H 0 ) 4 

-(H 2 \ (H 2 ) 4 
{H x \ -(H 3 ) 4 



(10) 



10 Le matrici (HO4, i = 0,1,2,3 hanno le seguenti espressioni: 



("0)4 = 



1000 
0100 
0010 
0001 



> (^1)4 = 



0 
j_ 

2 
0 



1 
0 

i_ 
2 
0 



0 

j_ 

2 
0 

j_ 
2 



0 
0 

1 
2 
0 
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0 0 1 



0 
2 



0 0 



0 
1_ 

2 



1 





0 0 



0 



0 



2 



La 1-D DCT e espressa dalla matrice (1-D DCT)4= 



1 



1 1 1 



C 4 = 




Q$ ~~Cs ~^8 

— c\ ~ c\ c\ 

~C£ ^8 ~^8 



, dove si e posto C™ = cos — 



Da queste equazioni si comprende che il calcolo di una DCT 4x4 
5 Pud essere suddiviso in due stadi: 

• il calcolo di quattro 1-D DCT, ognuna eseguita su una sequenza opportuna di 
quattro pixel 

• il calcolo della 2-D DCT a partire dalle quattro 1-D DCT 

Questi due stadi vengono eseguiti in modo simile, in modo tale da essere realizzati 
10 da un unico hardware, che viene utilizzato due volte. 

Si consideri ora come calcolare in parallelo quattro DCT 4x4. I 64 campioni 
complessivi sono ottenuti dai quattro blocchetti 4x4 in cui si suddivide un blocco 



II procedimento esuddiviso in varie fasi, a ciascuna della quali corrisponde un 
15 blocco architetturale. La veduta d'insieme e in figura 6. In questa figura vengono 
evidenziate le trasformazioni eseguite su ciascun blocchetto 4x4. 

Fase INPUT 

I pixel di ciascun quadrante (ij), 0 < /, j < 1 vengono ordinati per costituire i 



8x8. 
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vettori: 





■^0,0 




x l,0 




X 2,0 




x 3,0 


A \ - 


X l,l 

x 2,2 

x UJ 
_ 3,3 . 


, >»3 — 


x 3,l 
*0,2 
. 2.3 _ 


c'J _ 
> -°3 — 


x 0,l 
x 3,2 
_ x l,3 _ 


n'.y _ 
> a \ ~ 


JJ 
X 2,l 

x l,2 

x iJ 
. °.3 _ 



Dopo aver organizzato i dati in 16 vettori a 4 componenti, definiamo i vettori a 8 
componenti l,m,n t o 9 costituiti rispettivamente dalle componenti prima, seconda, 
terza, quarta degli otto vettori iniziali costituiti dai pixel dei quadranti 00 e 01 e i 
vettori p,q t r,s, costituiti rispettivamente dalle componenti prima, seconda, terza, 
quarta degli otto vettori iniziali costituiti dai pixel dei quadranti 10 e 11. 
Precisamente: 



A, [Of 0 " 




Adl]°'° 




Ad2T'° 
A 3 [ir 
B 3 [2f>° 




Ad3f fi 
A 3 [3f'° 
5 3 [3f° 
Bd3?'° 


43 [o] 0 - 0 

B 2 [0f>° 




A 3 [lf' 0 
B 2 [lf>° 














B X [if 0 




Bd2f>° 




Ad0f> 1 


, m = 


Adir 


, n = 


Ad*?' 1 


, o = 


Ad*?' 1 


A 3 [0f A 




A^r 








A 3 [3?' 1 


B 3 [of> 1 




B 3 [lf' 1 




B 3 [ 2 r 




B&r 


B l [Of 1 




M?' 1 . 




Bd2?' 1 







10 



p = 



Ad0f° 








~Ad2V~ 






A 3 [0V 




A 3 [lV 




A 3 [2r 




A 3 [3f'° 


B 3 [0] l '° 




bM>° 




B 3 [ 2 r 




5 3 [3p 


Bd0] l >° 
Ad0V 


, P = 


BdlV 
AdlV 




Bd2f° 
AdlV 




B&r 

Ad3V 


AM'' 




A 3 [ir 




A&r 
B&r 




A 3 [3]" 


BM' 1 




s 3 [ir 






B 3 [3V 


_Bd0f\ 








_Bd2f\ 




Mi] 1 ' 1 . 



Tenendo conto di come, a loro volta, sono stati definiti i vettori A[^ 9 A^ J ' 9 
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B l j J , B[' J si ottiene Pordiriamento mostrato in dettaglio in figura 7. Si osservi che 
in questa figura il blocco 8x8 originate e suddiviso nei quattro quadranti 4x4, 
nell'ambito di ciascun quadrante (ij), i pixel appartenenti ai vettori A^ ,J , A^ J , 

B*2 J , B[ iJ sono contrassegnati da colori diversi. 

5 Secondo quanto descritto sopra, il calcolo di una DCT 4x4 pud essere suddiviso in 
due stadi: di conseguenza, la fase PROCESS, che e Tunica fase in cui si svolgono 
operazioni aritmetiche, viene eseguita due volte: 

• la prima volta, per calcolare in parallelo le sedici 1-D DCT; 

• la seconda volta, per calcolare in parallelole quattro DCT 4x4 partendo dai 
10 coefficienti delle 1-D DCT, 

Per poter "sapere" se e in corso il primo o il secondo stadio del calcolo, si utilizza 
la variabile stage. 

Durante la fase INPUT la variabile stage viene aggiomata al valore 0. 

AlTingresso della fase PROCESS vi sono 64 MUX controllati dalla variabile 
15 stage. Ciascun MUX riceve due ingressi: 

• un pixel delPimmagine originate, proveniente dalla fase INPUT (questo 
ingresso viene selezionato quando stage — 0); 

• un coefficiente di una 1-D DCT, proveniente dalla fase ORDER (questo 
ingresso viene selezionato quando stage =1). 

20 Fase PROCESS 

Questa fase consiste nell'elaborare i vettori /, m, 5 come mostrato in figura 8. 

Nella figura 8 sono stati usati i seguenti simbolismi: 
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A = 



2C 8 x / 8x 3 

(#l) 4 
0 



0 

(#1)4 



per stage = 0 
per stage = 1 



B 



2C 8 j x/ gx8 

("3)4 
0 



0 

-(#3)4 



per stage = 0 
per stage = 1 



2Ci x /, 



8x8 



("2)4 



/?er stage = 0 
per stage = 1 



1 

2 



per .y/age = 0 
per stage = 1 5 



(11) 



(12) 



(13) 



(14) 



5 All'uscita della fase PROCESS vi sono 64 DEMUX controllati dalla variabile 
stage. I DEMUX hanno la funzione di indirizzare i dati secondo due possibilita: 

• se stage = 0 il dato in ingresso a ciascun DEMUX e un coefficiente di una 1-D 
DCT; pertanto il dato deve essere processato ulteriormente e, a tal fine, viene 
inviato alia fase ORDER; 
10 • se stage — 1 il dato in ingresso a ciascun DEMUX e un coefficiente di una 2-D 
DCT; pertanto il dato non deve essere processato ulteriormente e, a tal fine, 
viene inviato alia fase OUTPUT. 

Fase ORDER 

La fase ORDER consiste neU'ordinamento 



15 delle sequenze di uscita delle otto 1-D DCT in otto vettori m\.., s' cosi definiti: 
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Jo 

f 0,l 

Jo 



a[0] 
b[o] 
c[0] 
d[0] 
a[4] 
b[4] 
c[4] 
.4411 



f 0,0 

A 0 ' 1 





r«[2ji 








b[2] 








c[2] 






f 0,0 




d[2] 




70,0- 






a[6] 








b[6) 








c[6] 








_*1 











r -[031 










/[o] 
















" f l '°~ 
Jo 




h[0] 




\A l '°] 


Jo . 




e[4] 


_/i u 






/W 
*M 
_*[4L 







r = 





r*]i 






/[2] 






^[2] 




A' 0 




A[2] 










e[6] 








/[6] 






*[6] 
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Dopo la fase ORDER la variabile stage viene aggiomata al valore 1. I dati in 
uscita dalla fase ORDER vengono inviati alia fase PROCESS. 

Fase OUTPUT 

Questa fase consiste nel riordinamento dei dati provenienti dalla seconda (cioe con 
5 stage =1) esecuzione della fase PROCESS: a partire da questi dati, che 
costituiscono i vettori di 8 componenti a, b,.. t h y si costruisce il blocco 

di uscita Yn*n cosi definito: 



y[54] y[55] 



y[7l 

y[63\ 



-a[0] 


b[o] 


c[0] 


40] 


e\0] 


f[o] 


g[o] 


h[0]~ 


a[3] 


b[3] 


c[3] 


43] 


e[3] 


/[3] 


^[3] 


h[3] 


e\4] 


/W 


M 


*M 


a[4] 


b[4] 


c[4] 


d[4] 


_e[7] 


f[7] 


g[7] 


h[7] 


a[7] 


b[7] 


c[7] 


d[7\ 



(15) 



Le differenze principali fra l'hardware che realizza quattro DCT 4x4 e Thardware 
10 che calcola sedici DCT 2x2 sono le seguenti: 

• le sequenze in cui si devono ordinare i pixel di un blocco delPimmagine 
originale dipendono dalla dimensione della DCT richiesta; 

• per eseguire sedici DCT 2x2 la fase PROCESS deve essere eseguita una sola 
volta; invece per eseguire quattro DCT 4x4 la fase PROCESS deve essere 

1 5 eseguita due volte; 

• le operazioni eseguite nella fase PROCESS non sono sempre le stesse nei due 
casi. 

Calcolo di una DCT 8x8 

Per N = 8 TEq.(5) diventa: 



:SDOCID:<E1 9883052204> 



19 



J, J, r(2/+i>w "i r(2/+i> 

y>m n — / y X; / COS — 71 COS — ?T 



Ponendo: 



10 



dove: 





~yo~ 




7o" 


^64x1 = 


y\ 


> -^64x1 = 


/l 




yi. 




_/7. 



j'/'b'/.o yi,\ ■•■ ytjY 

{foM =Q =^r(k/to)' WL) = DC3r (K/to)' 

{^1,/ }J =0 = { x 0,0 x l,l x 2,2 x 3,3 x 4,4 x 5,5 x 6,6 x 7,7 

= ( x l,0 x 4,l x 7,2 x 5,3 x 2,4 x 0,5 x 3,6 x 6,7 

{ A S,i }J =0 = {*2,0 x 7,l x 3,2 x l,3 x 6,4 x 4,5 x 0,6 x 5,7 

{^7,/}J =0 = { x 3,0 x 5,l x l,2 x 7,3 x 0,4 x 6,5 x 2,6 x 4,7 

( S 7,i)if S5= o = x 2,l x 6,2 x 0,3 x 7,4 x l,5 x 5,6 x 3,7 
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fe,/}J =0 = {^5,0 x 0,l x 4,2 *6,3 x l,4 X 3,S x l,6 x 2,l}> 



fe,/ ij =0 = {*6,0 *3,1 *0,2 x 2,3 x 5,4 *7,5 x 4,6 x \j}> 

{ B lj} 7 i= 0 = {*7,0 x 6,l x 5,2 x 4 t 3 x 3,4 x 2,5 *1,6 x 0,l] 
si pud dimostrare che: 

*64xl = 0^)64 ^64x1 > 0 7 ) 



dove 



(* 8 )64 = 



(#o) 8 

(#l) 8 

("s) 8 



(# 3 ) 8 

-{H*X 
-(HlX 



(H 0 X 
(H 5 X 
-(H 6 X 
-(#l) 8 

-to) 8 

("?) 8 
(H 2 X 
(PlX 



(H Q X 

(HnX 

~iH 2 X 

~{HsX 
(H A X 
faX 

-fax 

~(HlX 



(HoX 

"(^ 7 ) 8 
-fe) 8 
(# 5 ) 8 
(# 4 ) 8 
-("3)8 
-(#6) 8 
("l) 8 



Le matrici (Hj)s> i=0,l,..,7 hanno le seguenti espressioni: 



("o) 8 = 



1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0 


1 



(#o) 8 

-(X5X 
(#l) 8 

-{HaX 

-to) 8 
fe) 8 



("o) 8 

-(" 3 ) 8 

freX 

{H n X 
-{n*X 

(PxX 
-{H2X 

(HsX 



{HqX 
-(#i) 8 
{h 2 X 

-{H 2 X 
(H4X 
~{HsX 

{H 6 X 
-{HnXi 



(18) 
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(#i) 8 = 



0 

l_ 
2 
0 

0 

0 

0 

0 

0 



0 

J_ 

2 
0 

0 

0 

0 

0 



0 0 0 0 

1 



0 0 0 0 

J_ 

2 
0 

J_ 

2 
0 



2 
0 

J_ 
2 

0 0 



0 
0 



0 
0 



0 0 

1 



2 
0 

l_ 
2 
0 



0 0 
1 



2 
0 

J_ 
2 
0 



(# 2 ) 8 = 



0 


0 


1 


0 


0 


0 


0 


0 


0 

1 

2 


1 

2 


0 


1 
2 


0 

1 

2 


0 


0 


0 


0 

1 

2 


0 


0 


0 

1 
2 


0 


0 


0 


0 

1 

2 


0 


0 


0 

1 

2 


0 


0 


0 


0 

1 

2 


0 


0 


0 

1 
2 


0 


0 


0 


0 

1 

2 


0 


0 


0 


0 


0 


0 


0 

1 

2 


0 


0 

1 

2. 


0 


0 


0 


0 


0 


0 
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lEtiiilciiiF 



(# 3 ) 8 = 



0 
0 

0 

J_ 
2 
0 

0 

0 

0 



0 

J_ 
2 
0 
J_ 
2 
0 



0 



0 

J_ 
2 
0 

0 

0 

J_ 
2 
0 

0 



1 

0 

0 

0 

0 

0 
J_ 
2 
0 



0 

1_ 
2 
0 

0 

0 

0 

0 

J_ 
2 



0 

0 

J_ 
2 
0 

0 

0 

0 

0 



0 
0 
0 

1_ 

2 
0 



0 
0 
0 

0 

J_ 
2 
0 



0 - 



2 



0 



0 
0 
0 
0 

j_ 

2 
0 

0 

0 



0 
J_ 
2 
0 

2 
0 

0 



0 

J_ 
2 
0 

0 

0 

1^ 

2 
0 



0 

2 
0 

0 

0 

0 

0 

2 



1 
0 
0 
0 
0 
0 
0 
0 



0 

J_ 

2 
0 

0 

0 

0 



0 
0 

J_ 

2 
0 

0 



2 



0 
0 
0 

j_ 

2 
0 



0 -- 



-- 0 



0 
0 
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0 
0 
0 
0 
0 

J_ 

2 
0 

0 



0 
0 



0 

j_ 

2 
0 

2 
0 



0 

l_ 
2 
0 



0 

J_ 
2 



0 

I 
2 
0 

0 

0 

0 

0 



0 

J_ 

2 
0 



0 
0 
0 



0 - 
2 



0 
0 
0 
0 



2 



I 
2 



0 
0 
0 

2 
0 

0 



0 
0 
0 
0 
0 
0 

1_ 

2 
0 



0 
0 
0 
0 
0 

J_ 

2 
0 



2 



0 
0 
0 

0 

J_ 
2 
0 

0 

0 



0 
0 

0 

]_ 
2 
0 

0 

0 



0 

I 

2 
0 

0 

0 



0 

J_ 
2 
0 

0 

0 



2 



-- 0 



0 



1 

0 
0 
0 

_i_ 

2 
0 

0 

0 
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{H 1 \ = 



0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 

l 
2 


i 
1 

2 


0 

-1 

1 

~2 


0 


0 


0 


0 


0 

1 

2 


0 

1 

2 


0 


0 


0 


0 

1 
2 


0 

_l 
2 


0 


0 


0 


0 

1 

2 


0 

1 
2 


0 


0 


0 


0 

1 

2 


0 

1 
2 


0 


0 


0 


0 

I 

2 


0 

_J_ 
2 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



La 1-D DCT e espressa dalla matrice 



(l-Z>DCr) 8 =C 8 = 



dove si e posto: 



1 


1 


1 


1 


1 


1 


1 


1 




ch 


cf. 


ch 


-ch 


-cf. 


-ch 


-ch 


c| 


c'i 


c| 


cl 


cl 


cl 


cl 


cl 




-ch 


-ch 


-ch 


ch 


ch 


ch 


-ch 


cj 


-cl 


-c\ 


ci 


ci 


-ci 


-ci 


cl 


cf. 


-cU 


ch 


CU 


-ch 


-ch 


ch 


-ch 


c| 


-cl 


cl 


ci 


ci 


cl 


cl 


cl 


.ch 


-ch 


ch 


-ch 


ch 


-ch 


-ch 


-ch. 



,m 



C™ =cos(— k) 
n 



Da queste equazioni si comprende che il calcolo di una DCT 8x8 pud essere 
suddiviso in due stadi: 



10 



• il calcolo di otto 1-D DCT, ognuna eseguita su una sequenza opportuna di otto 
pixel 

• il calcolo della 2-D DCT a partire dalle otto 1 -D DCT 



.HS^iS^PS 25 JB 
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Questi due stadi possono essere eseguiti in modo simile, in modo tale da essere 
realizzati da un unico hardware, che viene utilizzato due volte. 

II procedimento e suddiviso in varie fasi, a ciascuna della quali corrisponde un 
blocco architetturale. La veduta d'insieme e in figura 9. 

5 Fase INPUT 

I pixel del blocco 8x8 in ingresso vengono ordinati per costituire i vettori a 8 
componenti 1, m, n, o, p, q, r, s: 



>[0]' 








'A X m 


^ 3 [0] 




A 3 [l] 




A 3 m 


A 5 [0] 




A 5 [l] 




A 5 [71 


A 7 [0] 




A 7 [l] 




A 7 [7] 


Z? 7 [0] 


, m = 


B 7 [l] 




B 7 [7] 


B 5 [0] 




B 5 [l] 




B 5 [7] 


B 3 [0] 




B 3 [l] 




B 3 m 


5j[0] 











Tenendo conto di come, a loro volta, sono stati definiti i vettori Aj, A3, As, A 7, B7, 
10 B s , B 3 , Bj si ottiene Tordinamento mostrato in dettaglio in figura 10. Si osservi 
che in questa figura i pixel appartenenti ai vettori Aj t A3, As, Ay, By, Bs, B3, Bi 
sono contrassegnati da colori diversi. 

Come mostrato sopra, il calcolo di una DCT 8x8 pud essere suddiviso in due 
stadi: di conseguenza, la fase PROCESS, che e Tunica fase in cui si svolgono 
15 operazioni aritmetiche, viene eseguita due volte: 

• la prima volta, per calcolare in parallelo le sedici 1-D DCT; 

• la seconda volta, per calcolare la DCT 8x8 partendo dai coefficienti delle 1-D 
DCT. 

Per poter "sapere" se e in corso il primo o il secondo stadio del calcolo, si utilizza 
20 la variabile stage. 
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Durante la fase INPUT la variabile stage viene aggiornata al valore 0. 

All'ingresso della fase PROCESS vi sono 64 MUX controllati dalla variabile 
stage. Ciascun MUX riceve due ingressi: 

• un pixel deH'immagine originale, proveniente dalla fase INPUT (questo 
5 ingresso viene selezionato quando stage = 0); 

• un coefficiente di una 1-D DCT, proveniente dalla fase ORDER (questo 
ingresso viene selezionato quando stage = 1). 

Fase PROCESS 

Questa fase consiste neU'elaborare i vettori /, m, .... s come mostrato in figura 11. 
10 Nella figura 1 1 sono stati usati i seguenti simbolismi: 



A = 



B 



2Cg x / 8x8 per stage = 0 
{H 2 \ per stage = 1 



2Cg x / 8x8 per stage = 0 
-(H 6 \ per stage = 1 



(19) 



(20) 



c 2C\xI SxS per stage = 0 



(H 4 \ per stage = 1 



per stage = 0 
per stage = 1 



(22) 



15 All' uscita della fase PROCESS vi sono 64 DEMUX controllati dalla variabile 
stage. I DEMUX haimo la funzione di indirizzare i dati secondo due possibilita: 

• se stage = 0 il dato in ingresso a ciascun DEMUX e un coefficiente di una 1 -D 
DCT; pertanto il dato deve essere processato ulteriormente e, a tal fine, viene 
inviato alia fase ORDER; 
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• se stage = 1 il dato in ingresso a ciascun DEMUX e vm coefficiente della 2-D 
DCT; pertanto il dato non deve essere processato ulteriormente e, a tal fine, 
viene inviato alia fase OUTPUT. 



Fase ORDER 

La fase ORDER consiste nell'ordinamento delle sequenze di uscita delle ofto 1-D 
DCT in otto vettori m ', .... s' cosi definiti: 



am 
6[7] 
c[7] 

dm 

e[7] 
/[7] 

gm 

h[7] 



/'=/o = 



~a[0] 




"«[!]" 


b[0] 




b[l] 


c[0] 




c[l] 


d[0] 




d[l] 


e[0] 


e[l] 


/[0] 




/[I] 


£[0] 






>[0]_ 




_h[l] 



q' = f 7 = 



Dopo la fase ORDER la variabile stage viene aggiornata al valore 1. I dati in 
uscita dalla fase ORDER vengono inviati alia fase PROCESS. 

10 Fase OUTPUT 

Questa fase consiste nel riordinamento dei dati provenienti dalla seconda (cioe con 
stage =1) esecuzione della fase PROCESS: a partire da questi dati, che 
costituiscono i vettori di 8 componenti a f b, /r, si costruisce il blocco di uscita 
Yn+n cosi definito: 



15 



>>[54] y[55\ 



yVY 



a[0] Z>[0] c[0] d[0] e[0] /[0] g[0] MO] 

am bm cm dm em rm &m >m. 



Le differenze principali fra l'hardware che realizza una DCT 8x8 e l'hardware che 
calcola quattro DCT 4x4 so no le seguenti: 
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sequenze in cui si devono ordinare i pixel di un blocco dell'immagine 
originate dipendono dalla dimensione della DCT richiesta; 
• le operazioni eseguite nella fase PROCESS non sono sempre le stesse nei due 
casi. 

5 Procedimento di calcolo della DCT per blocchi di dimensione scalabile (DCT 
8x8, DCT 4x4 e DCT 2x1) 

Dai procedimenti gia esposti si ricava un algoritmo per calcolare, a scelta, una 
DCT 8x8 oppure 

quattro DCT 4x4 (in parallelo) oppure sedici DCT 2x2 (in parallelo). 

10 Questa scelta viene indicata dairutente assegnando il valore della variabile 
globale "size": 



size = 



0 per una DCT 8x8 

1 per quattro DCT 4x4 

2 per sedici DCT 2x2 



II procedimento e suddiviso in varie fasi (indipendentemente dal valore di "size "), 
a ciascuna della quali corrisponde un blocco architetturale. La veduta d'insieme e 
15 infigura 12. 

Ogni fase e stata riorganizzata in modo tale da fornire i risultati parziali 
corrispondenti al valore di prescelto, minimizzando le ridondanze. A volte le 
operazioni da eseguire sono diverse a seconda del valore di size: in questi casi 
rarchitettura prevede un MUX il cui ingresso di controllo e size. Si esaminino ora 
20 le varie fasi mettendo in evidenza le differenze rispetto alle architetture gia 
descritte: 

Fase INPUT 

La finalita di questa fase, descritta in Figura 13, 6 di ordinare i dati in modo 
opportuno per poter calcolare, partendo da loro, le 1-D DCT. Questo viene fatto 
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ricevendo in ingresso i valori di luminanza dei pixel (matrice x8x8), e ordinandoli 
negli otto vettori a otto componenti /, m, s. 

Per esempio: 



Fase PROCESS con staze = 0 

Questa fase consiste nel calcolare in parallelo le otto 1-D DCT elaborando i vettori 
/, ntj s come mostrato in figura 14. 

10 In questa figura si pud notare l'uso di 16 MUX controllati da size. 

Gli otto MUX sulla sinistra servono a poter "bypassare" le operazioni necessarie 
solo per il calcolo della DCT 8x8; pertanto il "bypass" awiene per size = 1 o 2, 
mentre non awiene per size = 0. 

Gli otto MUX sulla destra servono a inviare in uscita solo il risultato 
15 corrispondente al valore di size prescelto. 



l[0] = x 00 



5 




x 01 per size = 0 
s[7] = jx 4i7 per size = 1 
x 7 7 per size — 2 




2 



1 



per stage = 0 
per stage = 1 * 



(24) 
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2C g x / 8xg 



l 8x8 



0 

0 fa.V 



per stage, size = (0,0) oppure (0,1) 

per stage, size = (0,2) oppure (1,2) 
/?er stage, size = (1,0) 

/?er stage, size = (1,1) 



,(25) 



5 = 



0 



2C 8 x / 8x8 



^8x8 

0 

"(" 3 ) 4 



per stage, size — (0,0) oppure (0,1) 

/7er stage, size = (0,2) oppure (1,2) 
/?er stage, size = (1,0) 

/?er stage, size = (1,1) 



,(26) 



C = 



2Cj x /, 



8x8 



l 8x8 



(" 4 ) 8 



o (" 2 ) 4 



^?er stage, size = (0,0) oppure (0,1) 

/?er stage, size = (0,2) oppure (1,2) 
/?er stage, size = (1,0) 

/?er stage, size = (1,1) 



,(27) 



£> = 



2C'x/, 



16 ^ 8x8 



per 5tage = 0 
per stage = 1 



(28) 



5 E = 



2C 8 x I^g 



per stage = 0 
( H s\ per stage = \ 



(29) 



F = 



2ClxL 



16 ~ "* 8x8 



-(^ 7 ) 8 



/7er stage = 0 
per stage = 1 



(30) 



per stage = 0 



f 2C, 5 6 x/ 8)<8 

1 "C^Os per stage = 1 



(31) 



Lo schema in figura 14 pud essere scomposto nei blocchi mostrati in figura 15. 
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II blocco QA, per esempio, riceve in ingresso due vettori a otto componenti (ogni 
componente e un pixel, che pud aver gia subito elaborazioni) e invia come uscite 
due vettori a otto componenti: il primo vettore e la somma dei due vettori in 
ingresso, mentre il secondo vettore e la differenza fra i due vettori in ingresso 
5 successivamente eleborata tramite Toperatore lineare A. 

Si noti che gli operatori A,B,C,D,E,F,G sono matrici 8x8. 

Passando a un livello inferiore di generalizzazione i blocchi QA, QB, QC sono 
mo strati in dettaglio rispettivamente nelle figure 16, 17, 18. In queste figure i 
MUX sono controllati da tre bit, che corrispondono alia variabile stage (che pud 
10 valere 0 o 1, e pertanto e rappresentata da un bit) e alia variabile size (che puo 
valere 0,1 o 2, e pertanto e rappresentata da due bit). 

I blocchi QD, QE, QF, QG sono mostrati in dettaglio rispettivamente nelle figure 
19, 20, 21, 22. In queste figure i MUX sono controllati da un bit, che corrisponde 
alia variabile stage. 

15 Fase ORDER 

La fase ORDER, descritta in Figura 23, consiste nell'ordinamento delle sequenze 
di uscita delle otto 1-D DCT in otto vettori m\ s\ Per esempio: 



/'[0] = a[0]; 
/ , [1] = 6[0]; 



e[0] 



per size 



0 



20 



/'[4] = 



«[4] 
a[l] 



per size 
per size 



1 



2 
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47] = 47]; 



10 



Fase OUTPUT 

Questa fase, descritta in Figura 24, consiste nel riordinamento dei dati provenienti 
dalla seconda (cioe con stage = 1) esecuzione della fase PROCESS: a partire da 
questi dati, che costituiscono i vettori di 8 componenti a, b, ... , h si costruisce il 
blocco di uscita yNxN. 



Per esempio: 



M0] = 



a[0] per size — 0o\\ 

l[2] per size = 2 J 



y[l] = 



b[Q] per size = Ool] 

/[l] per size = 2 J 



y[63} = 



h[7] 
d[7] 
s[7] 



per size = 0 
per size = 1 
per size = 2 



Descrizione delle figure 

Uno schema funzionale a blocchi di un compressore-codificatore di immagini 
secondo la presente invenzione puo essere rappresentato come mostrato in Fizura 
15 7; Esso esegue pertanto una compressione "ibrida" che si basa essenzialmente su 
una codifica frattale eseguita nel dominio della DCT. Cio e reso possibile dalla 
peculiare architettura di calcolo parallelo della DCT su blocchi di pixel di 
dimensione scalabile, come descritto sopra. 

Vengono qui di seguito descritte singolarmente le altre figure. 

20 Fizura 2 : Grafo di flusso della DCT 2x2 
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Questo e il blocco "base', che compare ripetutamente nella fase PROCESS di 
tutte le DCT NxN, dove e una potenza di 2. 



In particolare: 

* II grafo di flusso per una DCT 2x2 e mostrato in figura 2, nella quale A = B = C 
5 = 1 e i dati in ingresso e in uscita sono i pixel nelle posizioni (0,0), (0,1 )", (1,0), 

* per sedici DCT 2x2 gli ingressi e le uscite sono vettori a 8 componenti e sono 
usati i seguenti simbolismi e si pone A = B = C = 18x8 

* per quattro DCT 4x4 gli ingressi e le uscite sono vettori a 8 componenti e sono 
10 usati i seguenti simbolismi: 



2Cl x /, 



S ~ * 8x8 



0 (#.)* 



per stage — 0 
per stage = 1 



(32) 



B = 



2C, 3 x/„ 



0 



0 



per stage = 0 
per stage = 1 



(33) 



C = 



2C\xI { 



8x8 



\H 2 \ 



per stage = 0 
per stage = 1 



(34) 



* per una DCT 8x8 gli ingressi e le uscite sono vettori a 8 componenti e sono usati 
15 i seguenti simbolismi: 



A = 



2Cg x 7 gx8 per stage = 0 



per stage = 1 



(35) 
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B = 



2C 8 x I SxS 



per stage — 0 
per stage - 1 



(36) 



C = 



2C\ x / 8x8 per stage = 0 



per stage = 1 



(37) 



* NelFarchitettura scalabile per il calcolo di una DCT 8x8 oppure quattro DCT 
4x4 (in parallelo) oppure sedici DCT 2x2 (in parallelo) gli ingressi e le uscite sono 
5 vettori a 8 componenti e sono usati i seguenti simbolismi: 



2C g x 



^8x8 



per (stage, size) = (0,0) oppure (0,1) 



per (stage, size) = (0,2) oppure (1,2) 



,(38) 



per (stage, size) (1,0) 



per (stage, size) (1,1) 



5 = 



2C 8 x i^g 



^ 8x8 
-(^) 8 

-(^3)4 0 
0 -(^ 3 ) 4 



per (stage, size) = (0,0) oppure (0,1) 



per (stage, size) = (0,2) oppure (1,2) 



/?er (stage, size) (1,0) 



/?er (stage, size) (1,1) 



(39) 
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iSPEG 



C = 



2C|x/ 8x8 



l 8x8 



0*4 X 

(^2)4 0 

L 0 (^2)4 



per (stage, size) = (0,0) oppure (0,1) 



per (stage, size) = (0,2) oppure (1,2) 



7?er (stage, size) (1,0) 



per (stage, size) (1,1) 



(40) 



Fizura 3 : Architettura per il calcolo di sedici DCT 2x2 

I pixel che costituiscono il blocco in ingresso vengono ordinati nella fase INPUT 
ed elaborati nella fase PROCESS in modo da ottenere i coefficient! delle sedici 2- 
5 D DCT su 4 campioni. Per esernpio, la 2-D DCT del blocchetto (0,1) costituito da 

{l[0] t m[0] t n[0],o[0]} e {a[0],b[0] M c[0Ld[0]}. 

I coefficienti delle 2-D DCT vengono riordinati nella fase ORDER in otto vettori a 
otto componenti. Per esernpio i coefficienti {a[0],b[0],c[0],d[0]} vanno a 
costituire il vettore / 

10 I sedici vettori a due componenti cosi ottenuti vengono inviati alia fase PROCESS 
per ottenere i coefficienti della DCT 2x2. Questi coefficienti, riordinati nella fase 
OUTPUT, costituiscono il blocco di uscita. 

Ft sura 4 : Ordinamento dei dati in ingresso per il calcolo di sedici DCT 2x2 

Questa figura mostra come i pixel del blocco 8x8 in ingresso vengono ordinati 
15 per costituire i vettori a 8 componenti 1, m, ... s. In ciascun quadrante (i,j), con 0 < 
i j < 3, i pixel appartenenti ai vettori sono contrassegnati da colori diversi. Per 
esernpio: 



Da ciascuno di questi vettori, le componenti con lo stesso indice (ossia i pixel con 
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mm 



lo stesso indice di colonna) vanno a formare uno stesso vettore a quattro 
componenti. Ad esempio il vettore 1 e costituito dagli elementi {A1[0], B1[0]}. 



Ne risulta che ciascun pixel del blocco 8x8 in ingresso viene a costituire una 
componente di uno dei vettori 1, m, n, o, p, q, r, s. 

5 Fizura 5 : Fase PROCESS per il calcolo di sedici DCT 2x2 

Questa fase consiste nell'elaborare i vettori a 8 componenti 1, m, , s. La fase 
PROCESS}, che e Tunica fase in cui si svolgono operazioni aritmetiche, viene 
eseguita una sola volta, per calcolare in parallelo le sedici 2-D DCT. 

Fizura 6 : Architettura per il calcolo di quattro DCT 4x4 

10 I pixel che costituiscono il blocco in ingresso vengono ordinati nella fase INPUT 
ed elaborati nella fase PROCESS in modo da ottenere i coefficienti delle sedici 1- 
D DCT su 4 campioni. Per esempio, la 1-D DCT della sequenza 
{l[Q],m[0],n[0],o[0]} i{a[0],b[0],c[0],df0]}. 

I coefficienti delle 1 -D DCT vengono riordinati nella fase ORDER in otto vettori a 
15 otto componenti. Per esempio i coefficienti {a[0],b[0],c[0],d[0]} vanno a 
costituire il vettore 1'. 

I quattro vettori a quattro componenti cosi ottenuti vengono inviati alia fase 
PROCESS per ottenere i coefficienti della DCT 4x4. Questi coefficienti, riordinati 
nella fase OUTPUT}, costituiscono il blocco di uscita. 

20 Figura 7 : Ordinamento dei dati in ingresso per il calcolo di quattro DCT 4x4 

Questa figura mostra come i pixel dei blocco 8x8 in ingresso vengono ordinati 
per costituire i vettori a 8 componenti 1, m, s. 

In ciascun quadrante (i j), con 0 < i j < 3, i pixel appartenenti ai vettori 
sono contrassegnati da colori diversi. Per esempio: 
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Da ciascuno di questi vettori, le componenti con lo stesso indice (ossia i pixel con 
lo stesso indice di colonna) vanno a formare uno stesso vettore a quattro 
componenti. Ad esempio il vettore 1 e costituito dagli elementi {A1[0], A3[0], 
5 B3[0], B1[0]}. 

Ne risulta che ciascun pixel del blocco 8x8 in ingresso viene a costituire una 
componente di uno dei vettori 1, m, n, o, p, q, r, s. 

Fizura 8 : Fase PROCESS per il calcolo di quattro DCT 4x4 

Questa fase consiste nelFelaborare i vettori a 8 componenti /, m, .... s. 

10 La fase PROCESS, che e Tunica fase in cui si svolgono operazioni aritmetiche, 
viene eseguita due volte: 

* la prima volta (stage = 0), per calcolare in parallelo le sedici 1-D DCT; 

* la seconda volta (stage = 1), per calcolare la DCT 8x8 partendo dai coefficienti 
delle 1-D DCT. 



15 A = 



2C 8 x 
("A 0 



per stage — 0 
per stage = 1 



(41) 



B = 



2C 8 x I^g 



0 



per stage = 0 
per stage — 1 



(42) 
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c = 



1C\ x I Sxi per stage = 0 

0 



(HA 0] f ' < 43 > 

per stage = 1 



(1 per stage = 0 

> (44) 
2 per stage = 1 

Fizura 9 : Architettura per il calcolo di una DCT 8x8 

I pixel che costituiscono il blocco in ingresso vengono ordinati nella fase INPUT 
5 ed elaborati nella fase PROCESS in modo da ottenere i coefficient! delle otto 1-D 
DCT su 8 campioni. Per esempio, la 1-D DCT della sequenza {l[0],m[0], ... ,s[0]} 
e{a[0],b[0], .... ,h[0]}. 

I coefficienti delle 1-D DCT vengono riordinati nella fase ORDER in otto vettori a 
otto componenti. Per esempio i coefficienti {a[0],b[0], .... ,h[7]} vanno a costituire 
10 il vettore 1'. 

Gli otto vettori a otto componenti cosi ottenuti vengono inviati alia fase 
PROCESS per ottenere i coefficienti della DCT 8x8. Questi coefficienti, riordinati 
nella fase OUTPUT, costituiscono il blocco di uscita. 

Figura 10 : Ordinamento dei dati in ingresso per il calcolo di una DCT 8x8 

15 Questa figura mostra come i pixel del blocco 8x8 in ingresso vengono ordinati 
per costituire i vettori a 8 componenti 1, m, .... , s. I pixel appartenenti ai vettori 
Al, A3, A5, A7 t B7 t B5, B3 t Bl sono contrassegnati da colori diversi. Per esempio: 



Da ciascuno di questi vettori, le componenti con lo stesso indice (ossia i pixel con 
20 lo stesso indice di colonna) vanno a formare uno stesso vettore a otto componenti. 
Ad esempio il vettore 1 e costituito dagli elementi { A1[0], A3[0], , B1[0]}. 
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Ne risulta che ciascun pixel del blocco 8x8 in ingresso viene a costituire una 
componente di uno dei vettori 1, m, n, o, p, q, r, s. 

Fizura 11 : Fase PROCESS per il calcolo di una DCT 8x8 

Questa fase consiste neirelaborare i vettori a 8 componenti 1, m, s. 

5 La fase PROCESS, che e Tunica fase in cui si svolgono operazioni aritmetiche, 
viene eseguita due volte: 

* la prima volta (stage — 0), per calcolare in parallelo le sedici 1-D DCT; 

* la seconda volta (stage = 1), per calcolare la DCT 8x8 partendo dai coefficienti 
delle 1-D DCT. 

10 Nella figura 1 1 sono stati usati i seguenti simbolismi: 

f 2C\*I per stage -0^ (45) 

I \H 2 ) S per stage = 1 

f 2C^x/ per stage = 0^ (46) 

I -{H 6 h per stage = I 

f 2C\xI per stage = 0 (4?) 

[ (tfJg per stage = I 

[ 1 per stage = ^ (4g) 

[ 2 per stage — 1 

15 Fizura 12 : Architettura scalabile per il calcolo di una DCT 8x8 o di quattro DCT 
4x4 o di sedici DCT 2x2 

I pixel che costituiscono il blocco in ingresso vengono ordinati nella fase INPUT 
ed elaborati nella fase PROCESS, dove vengono calcolate: 

* le DCT 1-D (per stage = 0, cioe per la DCT 8x8, e per stage = 1, cioe per le DCT 
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4x4) 

* direttamente le DCT 2-D (per stage} = 2, cioe per le DCT 2x2) 

Nei casi stage = 0e stage - 1 i coefficienti vengono poi riordinati nella fase 
ORDER in otto vettori a otto componenti, che vengono inviati alia fase PROCESS 
5 per ottenere i coefficienti della DCT 2-D. Questi coefficienti, riordinati nefla fase 
OUTPUT, costituiscono il blocco di uscita. 

Nel caso stage = 2 i coefficienti vengono trasmessi direttamente alia fase 
OUTPUT, dove vengono riordinati per costituire il blocco di uscita. 

Fizura 13 : Fase INPUT per l'architettura scalabile 

10 Gli ingressi sono i 64 pixel che costituiscono il blocco in ingresso. 

L'ordinamento degli ingressi viene operato per mezzo di MUX controllati dalla 
variabile size. 

Le 64 uscite sono le componenti degli otto vettori a otto componenti /, m, .... , 5. 

Fieura 14 : Fase PROCESS per l'architettura scalabile. 

15 Questa fase consiste nel calcolare in parallelo le otto 1-D DCT elaborando i vettori 
/, m, ... s come mostrato in figura 1 1 . 

In questa figura si pud notare Tuso di 16 MUX controllati da size. 

Gli otto MUX sulla sinistra servono a poter "bypassare' le operazioni necessarie 
solo per il calcolo della DCT 8x8; pertanto il "bypass» awiene per stage = 1 o2, 
20 mentre non awiene per stage = 0. 

Gli otto MUX sulla destra servono a inviare in uscita solo il risultato 
corrispondente al valore di size prescelto. 
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Nella Figura 14 sono stati usati i seguenti simbolismi: 



t = 



1 per stage = 0 

2 per stage = 1 3 



(49) 



2C 8 x 7 8x8 



l 8x8 



0 (»,) 4 



per (stage, size) = (0,0) oppure (0,1) - 



per {stage, size) = (0,2) oppure (1,2) 



,(50) 



per {stage, size) (1,0) 



per (stage, size) (1,1) 



2C 8 x 7 8x8 



^8x8 

-(^3)4 0 

0 -to). 



per (stage, size) = (0,0) oppure (0,1) 



per (stage, size) = (0,2) oppure (1,2) 



,(51) 



per (stage, size) (1,0) 



per (stage, size) (1,1) 



2C;x/ j 



8x8 



[ 8x8 



0 (tf 2 ) 4 



per (stage, size) = (0,0) oppure (0,1) 



per (stage, size) = (0,2) oppure (1,2) 



,(52) 



per {stage, size) (1,0) 



/?er {stage, size) (1,1) 



£> = ( 2C }* x I ** 



per stage = 0 
per stage = 1 



(53) 
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2Cj 3 6 x / M per stage = 0 
i H s\ per stage = V 



(54) 



[ -(H 7 \ per stage = T 

f 2C,W jurrf^-O^ (56) 
I -(^3)5 per stage = 1 

Figura 15 : Schema dei blocchi che implementano la fase PROCESS 

5 Blocchi che compongono la fase PROCESS per Tarchitettura scalabile. II blocco 
QA, per esempio, riceve in ingresso due vettori a otto componenti (ogni 
componente e un pixel, che pud aver gia subito elaborazioni) e invia come uscite 
due vettori a otto componenti: il primo vettore e la somma dei due vettori in 
ingresso, mentre il secondo vettore e la differenza fra i due vettori in ingresso 

10 successivamente elaborata tramite l'operatore lineare A. Si noti che gli operatori 
A,B,C,D,E,F,G sono matrici 8x8. 

Fiswra 16 : Schema dettagliato del blocco QA 

In questo schema sono mostrate in dettaglio le singole componenti dei due vettori 
in ingresso e gli operatori aritmetici (sommatori ecc.) che agiscono su ciascuna 
15 componente. I risultati vengono inviati ai MUX (sulla destra), ciascuno dei quali, 
in base ai valori delle variabili di controllo stage e size, seleziona un solo risultato, 
che costituisce una componente del vettore in uscita. 

Fizura 1 7 : Schema dettagliato del blocco QB 

In questo schema sono mostrate in dettaglio le singole componenti dei due vettori 
20 in ingresso e gli operatori aritmetici (sommatori ecc.) che agiscono su ciascuna 
componente. I risultati vengono inviati ai MUX (sulla destra), ciascuno dei quali, 
in base ai valori delle variabili di controllo stage 
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e size, seleziona un solo risultato, che costituisce vma componente del vettore in 
uscita. 

Figura 18 : Schema dettagliato del blocco QC 

In questo schema sono mostrate in dettaglio le singole componenti dei due vettori 
5 in ingresso e gli operatori aritmetici (sommatori ecc.) che agiscono su ciascuna 
componente. I risultati vengono inviati ai MUX (sulla destra), ciascuno dei quali, 
in base ai valori delle variabili di controllo stage 

e size, seleziona un solo risultato, che costituisce una componente del vettore in 
uscita. 

10 Figura 19 : Schema dettagliato del blocco QD 

In questo schema sono mostrate in dettaglio le singole componenti dei due vettori 
in ingresso e gli operatori aritmetici (sommatori ecc.) che agiscono su ciascuna 
componente. I risultati vengono inviati ai MUX (sulla destra), ciascuno dei quali, 
in base al valore della variabile di controllo stage, seleziona un solo risultato, che 
15 costituisce una componente del vettore in uscita. 

Figura 20 : Schema dettagliato del blocco QE 

In questo schema sono mostrate in dettaglio le singole componenti dei due vettori 
in ingresso e gli operatori aritmetici (sommatori ecc.) che agiscono su ciascuna 
componente. I risultati vengono inviati ai MUX (sulla destra), ciascuno dei quali, 
20 in base al valore della variabile di controllo stage, seleziona un solo risultato, che 
costituisce una componente del vettore in uscita. 

Figura 21 : Schema dettagliato del blocco QF. In questo schema sono mostrate in 
dettaglio le singole componenti dei due vettori in ingresso e gli operatori 
aritmetici (sommatori ecc.) che agiscono su ciascuna componente. I risultati 
25 vengono inviati ai MUX (sulla destra), ciascuno dei quali, in base al valore della 
variabile di controllo {\em stage}, seleziona un solo risultato, che costituisce una 
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componente del vettore in uscita. 

Fizura 22 : Schema dettagliato del blocco QG. In questo schema sono mostrate in 
dettaglio le singole cornponenti dei due vettori in ingresso e gli operatori 
aritmetici(sommatori ecc.) che agiscono su ciascuna componente. I risultati 
5 vengono inviati ai MUX (sulla destra), ciascuno dei quali, in base al valore della 
variabile di controllo stage, seleziona un solo risultato, che costituisce una 
componente del vettore in uscita. 

Fizura 23 : Fase ORDER per Parchitettura scalabile. Gli ingressi sono i 64 pixel 
dopo essere stati elaborati nella fase PROCESS. 

10 L'ordinamento degli ingressi viene operato per mezzo di MUX controllati dalla 
variabile size. 

Le 64 uscite sono le cornponenti degli otto vettori a otto cornponenti I', m\ s\ 

Fisitra 24 : Fase OUTPUTper Parchitettura scalabile.Gli ingressi sono i 64 
coefficient! 2-D DCT. 

15 L'ordinamento degli ingressi viene operato per mezzo di MUX controllati dalla 
variabile size. 

Le 64 uscite sono i pixel che costituiscono il blocco di uscita. 
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RIVENDICAZIONI 



1 . Metodo di calcolo della trasformata coseno discreta (DCT) per blocchi 
di pixels di un'immagine, caratterizzato dal fatto che comprende le operazioni di 
definire primi blocchi di suddivisione di detta immagine denominati blocchi 
5 range, aventi una dimensione frazionaria e scalabile di N/2'*N/2 l , dove"/ e un 
numero intero, rispetto ad una dimensione massima predefinita di un blocco di 
N*N pixels, denominato blocco domain, traslabile per intervalli di N/2* pixels, ed 
effettuare il calcolo della DCT parallelamente su 2' blocchi range di suddivisione 
di un blocco domain di N*N pixels. 

10 2. Metodo secondo la rivendicazione 1, caratterizzato dal fatto che il 

calcolo della DCT parallelamente su tutti i blocchi range di suddivisione di un 
certo blocco domain e eseguito mediante una struttura hardware e comprende le 
seguenti operazioni: 

a) ordinare i pixels in funzione della suddivisione in blocchi range di una certa 
15 dimensione riordimando i pixels di ingresso in un numero 2' di sequenze o 

vettori di 2' componenti; 

b) calcolare in parallelo 2 l DCT monodimensionali elaborando detti vettori 
definiti nella precedente fase a); 

c) ordinare le sequenze di uscita dei dati monodimensionali DCT relativi a 
20 detti 2' vettori; 

d) completare il calcolo parallelo di detti 2 l DCT bidimensionali elaborando 
dette sequenze di uscita prodotte nella fase c); 

e) ordinare le sequenze di uscita prodotte nella fase d) in un numero 2' di 
vettori di coefficienti DCT bidimensionali. 

25 3. Metodo secondo la rivendicazione 2, caratterizzato dal fatto che il 

calcolo parallelo di dette 2' DCT monodimensionali della fase b) e il 
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completamento del calcolo parallelo di dette 2 l DCT tridimensional! della fase d) 
sono effettuati suddividendo le sequenze risultanti rispettivamente dalla fase a) e 
dalla fase c) in gruppi di quattro elementi scalari, calcolandone le relative somme 
e differenze mediante sommatori e sottrattori e moltiplicando i risultati somma e 
5 differenza per rispettivi coefficienti in modo reiterativo fino a cornpletare il 
calcolo dei relativi quattro coefficienti DCT, rispettivamente monodimensionali e 
Tridimensional!. 

4. Metodo di compressione dei dati di un'immagine, da memorizzare o 
trasmettere, mediante codifica frattale, caratterizzato dal fatto che la 
trasfonnazione frattale e effettuata nel domain della trasformata coseno discreta 
(DCT) attraverso le seguenti operazioni 

• suddividere un'immagine in detti due distinti tipi di blocchi di pixels; 

• calcolare parallelamente la trasformata coseno discreta (DCT) di tutti i 2 l 
blocchi range e di un relativo blocco domain; 

• classificare i blocchi range trasformati in funzione della loro relativa 
complessita rappresentata dalla somma dei valori della terna di coefficienti AC; 

• applicare la trasformata frattale nel dominio della DCT ai dati dei blocchi 
range con classificazione di complessita superiore ad una soglia prestabilita e 
memorizzare solo il coefficiente DC dei blocchi range con complessita 
inferiore a detta soglia, identificando un relativo blocco domain di 
appartenenza del blocco range in trasfonnazione tale da produrre la migliore 
approssimazione frattale dello stesso blocco range; 

• calcolare un'immagine differenza tra ciascun blocco range e la sua 
approssimazione frattale; 

25 • quantizzare detta immagine differenza nel dominio della DCT impiegando una 
tabella di quantizzazione predisposta in funzione delle caratteristiche della vista 
umana; 

i SDOCID:<E1 98S30522O4> 





• codificare rimmagine differenza quantizzata mediante un metodo di codifica 
basato sulla probability dei coefficienti di quantizzazione; 

• memorizzare o trasmettere il codice di codifica per ciascun blocco range 
compresso nel dominio della DCT e il coefficiente DC di ciascun blocco range 

5 non compresso. 
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